Variant- and Vaccination-Specific Alternative Splicing Profiles in SARS-CoV-2 Infections (preprint)

The COVID-19 pandemic, caused by the coronavirus SARS-CoV-2, and its subsequent variants has underscored the importance of understanding the host-viral molecular interactions to devise effective therapeutic strategies. A significant aspect of these interactions is the role of alternative splicing in modulating host responses and viral replication mechanisms. Our study sought to delineate the patterns of alternative splicing of RNAs from immune cells across different SARS-CoV-2 variants and vaccination statuses, utilizing a robust dataset of 190 RNA-seq samples from our previous studies, encompassing an average of 212 million reads per sample. We identified a dynamic alteration in alternative splicing and genes related to RNA splicing were highly deactivated in COVID-19 patients and showed variant- and vaccination-specific expression profiles. Overall, Omicron-infected patients exhibited a gene expression profile akin to healthy controls, unlike the Alpha or Beta variants. However, significantly, we found identified a subset of infected individuals, most pronounced in vaccinated patients infected with Omicron variant, that exhibited a specific dynamic in their alternative splicing patterns that was not widely shared amongst the other groups. Our findings underscore the complex interplay between SARS-CoV-2 variants, vaccination-induced immune responses, and alternative splicing, emphasizing the necessity for further investigations into these molecular cross-talks to foster deeper understanding and guide strategic therapeutic development.

Blood transcriptomics analysis offers insights into variant-specific immune response to SARS-CoV-2 (preprint)

Bulk RNA sequencing (RNA-seq) of blood is typically used for gene expression analysis in biomedical research but is still rarely used in clinical practice. In this study, we argue that RNA-seq should be considered a routine diagnostic tool, as it offers not only insights into aberrant gene expression and splicing but also delivers additional readouts on immune cell type composition as well as B-cell and T-cell receptor (BCR/TCR) repertoires. We demonstrate that RNA-seq offers vital insights into a patients immune status via integrative analysis of RNA-seq data from patients infected with various SARS-CoV-2 variants (in total 240 samples with up to 200 million reads sequencing depth). We compare the results of computational cell-type deconvolution methods (e.g., MCP-counter, xCell, EPIC, quanTIseq) to complete blood count data, the current gold standard in clinical practice. We observe varying levels of lymphocyte depletion and significant differences in neutrophil levels between SARS-CoV-2 variants. Additionally, we identify B and T cell receptor (BCR/TCR) sequences using the tools MiXCR and TRUST4 to show that - combined with sequence alignments and pBLAST - they could be used to classify a patients disease. Finally, we investigated the sequencing depth required for such analyses and concluded that 10 million reads per sample is sufficient. In conclusion, our study reveals that computational cell-type deconvolution and BCR/TCR methods using bulk RNA-seq analyses can supplement missing CBC data and offer insights into immune responses, disease severity, and pathogen-specific immunity, all achievable with a sequencing depth of 10 million reads per sample. Key PointsO_LIComputational deconvolution of transcriptomes can estimate immune cell abundances in SARS-CoV-2 patients, supplementing missing CBC data. C_LIO_LI10 million RNA sequencing reads per sample suffice for analyzing immune responses and disease severity, including BCR/TCR identification. C_LI

Comprehensive analysis of immune responses in CLL patients after heterologous COVID-19 vaccination (preprint)

Patients with chronic lymphocytic leukemia (CLL) treated with B-cell pathway inhibitors and anti-CD20 antibodies exhibit low humoral response rate (RR) following SARS-CoV-2 vaccination. To investigate the relationship between the initial transcriptional response to vaccination with ensuing B and T cell immune responses, we performed a comprehensive immune transcriptome analysis flanked by antibody and T cell assays in peripheral blood prospectively collected from 15 CLL/SLL patients vaccinated with heterologous BNT162b2/ChAdOx1 with follow up at a single institution. The two-dose antibody RR was 40% increasing to 53% after booster. Patients on BTKi, venetoclax ± anti-CD20 antibody within 12 months of vaccination responded less well than those under BTKi alone. The two-dose T cell RR was 80% increasing to 93% after booster. Transcriptome studies revealed that seven patients showed interferon-mediated signaling activation within 2 days and one at 7 days after vaccination. Increasing counts of COVID-19 specific IGHV genes correlated with B-cell reconstitution and improved humoral RR. T cell responses in CLL patients appeared after vaccination regardless of treatment status. A higher humoral RR was associated with BTKi treatment and B-cell reconstitution. Boosting was particularly effective when intrinsic immune status was improved by CLL-treatment.

Omicron BA.1 and BA.2 immune response in naive and prior infected persons (preprint)

The highly transmissible SARS-CoV-2 Omicron (B.1.1.529) variant has replaced previous variants and is less susceptible to neutralizing antibodies elicited by vaccination or infection. Currently, BA.2 is the dominant omicron sublineage. Vaccinated individuals with BA.1 infection develop measurable neutralizing antibody titers against BA.1 and BA.2. The ability of BA.2 infection to induce neutralizing antibodies in unvaccinated individuals, either without or with previous SARS-CoV-2 infection, is pending definition.

Prior vaccination enables a more robust immune response to Omicron infection (preprint)

The antibody response following Omicron infection has been reported to be less robust than to other variants. Here we compared the immune-transcriptome and antibody responses following Omicron infection in unvaccinated and vaccinated individuals who experienced mild to moderate symptoms. The unvaccinated individuals showed a quantitatively greater transcriptional response but a muted antibody response than vaccinated individuals. Prior vaccination modifies the transcriptional response to Omicron infection with a more robust antibody response.

Genetic immune response and antibody repertoire of heterologous ChAdOx1-BNT162b2 vaccination in a Korean cohort (preprint)

Heterologous ChAdOx1-BNT162b2 vaccination induces a stronger immune response than two doses of BNT162b2 or ChAdOx1. Yet, the molecular transcriptome, the germline allelic variants of immunoglobulin loci and anti-Omicron antibody levels induced by the heterologous vaccination have not been formally investigated. Moreover, there is a paucity of COVID vaccine studies including diverse genetic populations. Here, we show a robust molecular immune transcriptome and antibody repertoire in 51 office workers from the Republic of Korea after a heterologous ChAdOx1-BNT162b2 vaccination or a homologous ChAdOx1-ChAdOx1 vaccination. Anti-spike-specific IgG antibody levels in the heterologous group increased from 14,000 U/ml to 142,000 AU/ml within eight days after the BNT162b2 vaccination. In contrast, antibody levels in the homologous group increased two-fold after the second ChAdOx1 dose. Antibody titers against the Omicron spike protein as compared to the ancestral strain were reduced to a lesser extent in the heterologous group. RNA-seq conducted on immune cells demonstrated a stronger activation of interferon-induced genetic programs in the heterologous cohort. An increase of specific IGHV clonal transcripts encoding neutralizing antibodies was preferentially detected in the heterologous cohort. Enrichment of B cell and CD4+ T cell responses were observed following both heterologous and homologous vaccination using scRNA-seq, but clonally expanded memory B cells were relatively stronger in the ChAdOx1-BNT162b2 cohort. In summary, a heterologous vaccination with ChAdOx1 followed by BNT162b2 provides an innate and adaptive immune response exceeding that seen in homologous ChAdOx1 vaccinations but equivalent to that seen in homologous BNT162b2 vaccination.

Robust immune response to the BNT162b mRNA vaccine in an elderly population vaccinated 15 months after recovery from COVID-19 (preprint)

Knowledge about the impact of prior SARS-CoV-2 infection of the elderly on mRNA vaccination response is needed to appropriately address the need for booster vaccination in this vulnerable population. To address this, we investigated antibody and genomic immune responses in 16 elderly (avg. 81 yrs.) individuals that had received a single booster dose of BNT162b vaccine 15 months after recovering from COVID-19. Spike-specific IgG antibody levels increased in each of the study participants from an average of 710 U/ml prior to the vaccination to more than 40,000 U/ml within ten weeks after the vaccination. In contrast, anti-spike-specific IgG antibody levels averaged 2,190 U/ml in 14 healthy SARS-CoV-2-naive individuals (avg. 58 yrs.) ten weeks after the second dose of BNT162b. RNA-seq conducted on PBMCs demonstrated the activation of interferon-activated genetic programs in both cohorts within one day. Unlike their transient induction in the younger naive population, persistent activity and the initiation of additional cell cycle regulated programs were obtained in the older COVID-19 recovered population. Here we show that the elderly, a high-risk population, can mount a strong antibody and a persistent molecular immune response upon receiving a single dose of mRNA vaccine 15 months after recovery from COVID-19.

Immune transcriptomes from hospitalized patients infected with the SARS-CoV-2 variants B.1.1.7 and B.1.1.7 carrying the E484K escape mutation (preprint)

Fast-spreading variants of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) energize the COVID-19 pandemic. B.1.1.7 (VOC-202012/01) has become the predominant variant in many countries and a new lineage (VOC-202102/02) harboring the E484K escape mutation in the B.1.1.7 background emerged in February 2021. This variant is of concern due to reduced neutralizing activity by vaccine-elicited antibodies. However, it is not known whether this single amino acid change leads to an altered immune response. Here, we investigate differences in the immune transcriptome in hospitalized patients infected with either B.1.1.7 (n=28) or B.1.1.7+E484K (n=12). RNA-seq conducted on PBMCs isolated within five days after the onset of COVID symptoms demonstrated elevated activation of specific immune pathways, including JAK-STAT signaling, in B.1.1.7+E484K patients as compared to B.1.1.7. Longitudinal transcriptome studies demonstrated a delayed dampening of interferon-activated pathways in B.1.1.7+E484K patients. Prior vaccination with BNT162b vaccine (n=8 one dose; n=1 two doses) reduced the transcriptome inflammatory response to B.1.1.7+E484K infection relative to unvaccinated patients. Lastly, the immune transcriptome of patients infected with additional variants (B.1.258, B.1.1.163 and B.1.7.7) displayed a reduced activation compared to patients infected with B.1.1.7. Acquisition of the E484K substitution in the B.1.1.7 background elicits an altered immune response, which could impact disease progression.

Impact of BNT162b first vaccination on the immune transcriptome of elderly patients infected with the B.1.351 SARS-CoV-2 variant (preprint)

Fast-spreading variants of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) energize the COVID-19 pandemic. The B.1.351 variant carrying the escape mutation E484K in the receptor binding domain is of particular concern due to reduced immunological protection following vaccination. Protection can manifest as early as 10 days following immunization with full protection two weeks following the second dose, but the course is not well-characterized for variants. Here, we investigated the immune transcriptome of six elderly individuals (average age 82 yr.) from an old peoples home, who contracted B.1.351, with four having received the first dose of BNT162b eight to 11 days prior to the onset of COVID-19 symptoms. The patients were hospitalized and received dexamethasone treatment. Immune transcriptomes were established from PBMCs approximately 10 and 35 days after the onset of COVID-19 symptomology. RNA-seq revealed a more intensive immune response in vaccinated patients as compared to unvaccinated ones. Specifically, transcription factors linked to the JAK/STAT pathway, interferon stimulated genes, and genes associated with innate antiviral immunity and COVID-19-SARS-CoV-2 infection were highly enriched in vaccinated patients. This rendered the transcriptomes of the older vaccinated group significantly different than older unvaccinated individuals infected at the same institution and more similar to the immune response of younger unvaccinated individuals (age range 48-62) following B.1.351 infection. All individuals in this study whether vaccinated or not were hospitalized due to B.1.351 infection and one vaccinated patient died illustrating that although an enhanced immune response was documented infection it was insufficient to protect from disease. This highlights the need for maintaining physical distancing and prevention measures throughout the time course of vaccination in older adults.

Impact of BNT162b First Vaccination on the Immune Transcriptome of Elderly Patients Infected with the B.1.351 SARS-CoV-2 Variant (preprint)

Fast-spreading variants of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) energize the COVID-19 pandemic. The B.1.351 variant carrying the escape mutation E484K in the receptor binding domain is of particular concern due to reduced immunological protection following vaccination. Protection can manifest as early as 10 days following immunization with full protection two weeks following the second dose, but the course is not well-characterized for variants. Here, we investigated the immune transcriptome of six elderly individuals (average age 82 yr.) from an old people’s home, who contracted B.1.351, with four having received the first dose of BNT162b eight to 11 days prior to the onset of COVID-19 symptoms. The patients were hospitalized and received dexamethasone treatment. Immune transcriptomes were established from PBMCs approximately 10 and 35 days after the onset of COVID-19 symptomology. RNA-seq revealed a more intensive immune response in vaccinated patients as compared to unvaccinated ones. Specifically, transcription factors linked to the JAK/STAT pathway, interferon stimulated genes, and genes associated with innate antiviral immunity and COVID-19-SARS-CoV-2 infection were highly enriched in vaccinated patients. This rendered the transcriptomes of the older vaccinated group significantly different than older unvaccinated individuals infected at the same institution and more similar to the immune response of younger unvaccinated individuals (ages 48-62) following B.1.351 infection. All individuals in this study whether vaccinated or not were hospitalized due to B.1.351 infection and one vaccinated patient died illustrating that although an enhanced immune response was documented infection it was insufficient to protect from disease. This highlights the need for maintaining physical distancing and prevention measures throughout the time course of vaccination in older adults.

Activation of Interferon-Stimulated Transcriptomes and ACE2 Isoforms in Human Airway Epithelium Is Curbed by Janus Kinase Inhibitors (preprint)

SARS-CoV-2 infection of human airway epithelium activates genetic programs that lead to progressive hyperinflammation in COVID-19 patients. Here we report on genetic programs activated by interferons and the suppression by Janus kinase (JAK) inhibitors. The angiotensin-converting enzyme 2 (ACE2) is the receptor for SARS-CoV-2 and deciphering its regulation is paramount for understanding the cell tropism of SARS-CoV-2 infection. We identified candidate regulatory elements in the ACE2 locus in human primary airway cells and lung tissue. Activating histone and promoter marks and Pol II loading characterize the intronic dACE2 and define novel candidate enhancers distal to the genuine ACE2 promoter and within additional introns. dACE2, and to a lesser extent ACE2, RNA levels increased in primary cells treated with interferons and this induction was mitigated by JAK inhibitors that are used therapeutically in COVID-19 patients. Our analyses provide insight into ACE2 regulatory elements and highlight JAK inhibitors as suitable tools to suppress interferon-activated genetic programs in bronchial cells.

Activation of ACE2 and interferon-stimulated transcriptomes in human airway epithelium is curbed by Janus Kinase inhibitors (preprint)

SARS-CoV-2 infection of human airway epithelium activates genetic programs that lead to progressive hyperinflammation in COVID-19 patients. Here we report on the transcriptomic response of airway epithelium to interferons and its suppression by the JAK inhibitors Baricitinib and Ruxolitinib. There is a debate on the regulation of the conventional versus the novel intronic promoter inducing the short ACE2 isoform. Through RNA-seq and ChIP-seq analyses for activating chromatin marks and Polymerase II, we define the interferon-activated intronic regulatory region. Our results also support that the conventional ACE2 promoter is controlled by interferon.

Prime-boost protein subunit vaccines against SARS-CoV-2 are highly immunogenic in mice and macaques (preprint)

SARS-CoV-2 vaccines are advancing into human clinical trials, with emphasis on eliciting high titres of neutralising antibodies against the viral spike (S). However, the merits of broadly targeting S versus focusing antibody onto the smaller receptor binding domain (RBD) are unclear. Here we assessed prototypic S and RBD subunit vaccines in homologous or heterologous prime-boost regimens in mice and non-human primates. We find S is highly immunogenic in mice, while the comparatively poor immunogenicity of RBD was associated with limiting germinal centre and T follicular helper cell activity. Boosting S-primed mice with either S or RBD significantly augmented neutralising titres, with RBD-focussing driving moderate improvement in serum neutralisation. In contrast, both S and RBD vaccines were comparably immunogenic in macaques, eliciting serological neutralising activity that generally exceed levels in convalescent humans. These studies confirm recombinant S proteins as promising vaccine candidates and highlight multiple pathways to achieving potent serological neutralisation.

Immune transcriptomes of highly exposed SARS-CoV-2 asymptomatic seropositive versus seronegative individuals from the Ischgl community (preprint)

To investigate prevalence of ongoing activation of inflammation following asymptomatic SARS-CoV-2 infection we characterized immune cell transcriptomes from 43 asymptomatic seropositive and 52 highly exposed seronegative individuals with few underlying health issues following a community superspreading event. Four mildly symptomatic seropositive individuals examined three weeks after infection as positive controls demonstrated immunological activation. Approximately four to six weeks following the event, the two asymptomatic groups showed no significant differences. Two seropositive patients with underlying genetic disease impacting immunological activation were included (Cystic Fibrosis (CF), Nuclear factor-kappa B Essential Modulator (NEMO) deficiency). CF, but not NEMO, associated with significant immune transcriptome differences including some associated with severe SARS-CoV-2 infection (IL1B, IL17A, respective receptors). All subjects remained in their usual state of health from event through five-month follow-up. Here, asymptomatic infection resolved without evidence of prolonged immunological activation. Inclusion of subjects with underlying genetic disease illustrated the pathophysiological importance of context on impact of immunological response.

Immune transcriptomes of highly exposed SARS-CoV-2 asymptomatic seropositive versus seronegative individuals from the Ischgl community (preprint)

To investigate prevalence of ongoing activation of inflammation following asymptomatic SARS-CoV-2 infection we characterized immune cell transcriptomes from 43 asymptomatic seropositive and 52 highly exposed seronegative individuals with few underlying health issues following a community superspreading event. Four mildly symptomatic seropositive individuals examined three weeks after infection as positive controls demonstrated immunological activation. Approximately four to six weeks following the event, the two asymptomatic groups showed no significant differences. Two seropositive patients with underlying genetic disease impacting immunological activation were included (Cystic Fibrosis (CF), Nuclear factor-kappa B Essential Modulator (NEMO) deficiency). CF, but not NEMO, associated with significant immune transcriptome differences including some associated with severe SARS-CoV-2 infection (IL1B, IL17A, respective receptors). All subjects remained in their usual state of health from event through five-month follow-up. Here, asymptomatic infection resolved without evidence of prolonged immunological activation. Inclusion of subjects with underlying genetic disease illustrated the pathophysiological importance of context on impact of immunological response.

Activation of the SARS-CoV-2 Receptor Ace2 by Cytokines Through Pan JAK-STAT Enhancers (preprint)

ACE2, in concert with the protease TMPRSS2, binds the novel coronavirus SARS-CoV-2 and facilitates its cellular entry. The ACE2 gene is expressed in SARS-CoV-2 target cells, including Type II Pneumocytes (Ziegler, 2020), and is activated by interferons. Viral RNA was also detected in breast milk (Wu et al., 2020), raising the possibility that ACE2 expression is under the control of cytokines through the JAK-STAT pathway. Here we show that Ace2 expression in mammary tissue is induced during pregnancy and lactation, which coincides with the establishment of a candidate enhancer. The prolactin-activated transcription factor STAT5 binds to tandem sites that coincide with activating histone enhancer marks and additional transcription components. The presence of pan JAK-STAT components in mammary alveolar cells and in Type II Pneumocytes combined with the autoregulation of both STAT1 and STAT5 suggests a prominent role of cytokine signaling pathways in cells targeted by SARS-CoV-2.Funding: This work was supported by the Intramural Research Program (IRP) of the National Institute of Diabetes and Digestive and Kidney Diseases (NIDDK)and utilized the computational resources of the NIH HPC Biowulf cluster (http://hpc.nih.gov)Declaration of Interest: The authors declare not competing interests.

Activation of the SARS-CoV-2 receptor Ace2 by cytokines through pan JAK-STAT enhancers (preprint)

ACE2, in concert with the protease TMPRSS2, binds the novel coronavirus SARS-CoV-2 and facilitates its cellular entry. The ACE2 gene is expressed in SARS-CoV-2 target cells, including Type II Pneumocytes (Ziegler, 2020), and is activated by interferons. Viral RNA was also detected in breast milk (Wu et al., 2020), raising the possibility that ACE2 expression is under the control of cytokines through the JAK-STAT pathway. Here we show that Ace2 expression in mammary tissue is induced during pregnancy and lactation, which coincides with the establishment of a candidate enhancer. The prolactin-activated transcription factor STAT5 binds to tandem sites that coincide with activating histone enhancer marks and additional transcription components. The presence of pan JAK-STAT components in mammary alveolar cells and in Type II Pneumocytes combined with the autoregulation of both STAT1 and STAT5 suggests a prominent role of cytokine signaling pathways in cells targeted by SARS-CoV-2.